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THREE-DIMENSIONAL FINITE-ELEMENT ELASTIC ANALYSIS OF A 


THERMALLY CYCLED SINGLE-EDGE WEDGE GEOMETRY SPECIMEN 

by Peter T. Bizon*, Richard J. Hill**, Bruce P. Guilliams*, 

Sandra K. Drake**, and Jeffrey L. Kladden** 

*National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 

**Air Force Aero Propulsion Laboratory 
Wright-Patterson Air Force Base, Ohio 

SUMMARY 

Conventional three-dimensional finite-element elastic analysis techniques 
were used to analyze a single-edge wedge geometry specimen thermally cycled 
in fluidized beds . This analysis was performed as a joint program of the Lewis Re- 
search Center of NASA and the Air Force Aero Propulsion Laboratory (AFAPL) . Lewis 
used a general structural analysis program (NASTRAN) with a fine mesh to analyze a 
severe heating time increment of one alloy . The AFAPL used a specialized blade and 
disk program (IS03DQ) and a coarse mesh for the same time increment and alloy . 
Comparison of results for this combination using both methods showed very good 
agreement . The remaining combinations of test conditions and alloys considered 
were then analyzed using the coarse mesh model with the IS03DQ computer program . 

Seven different combinations of alloys and cycling conditions were analyzed 
using specimens of identical geometry. The alloys were NASA TAZ-8A, 316 stainless 
steel, and A-286, The cycling condition for the first mentioned alloy was alternate 
3-minute immersions in beds maintained at 316° and 1088° C (600° and 1990° F) . 
Cycling conditions for each of the other two alloys were alternate 4-minute immer- 
sions in a cooling bed always maintained at 21° C (70° F) and a heating bed main- 
tained at 740°, 850°, or 960° C (1364°, 1562°, or 1760° F). The test conditions se- 
lected were those that would result in reasonable failure times for the alloys . The 
most severe testing condition was used for the high-temperature alloy (NASA TAZ-8A) 
with less severe conditions used for the other two intermediate-temperature alloys . 
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The temperature , strain and stress at the critical location for a typical cycle 
are presented herein for each of the seven combinations analyzed. In addition, the 
temperature , strain and stress distributions over the complete wedge are given at 
the time of both maximum and minimum longitudinal leading-edge strain . The maxi- 
mum longitudinal strainrange (algebraic difference between maximum and minimum 
longitudinal strain) varied from 0.35 to 0.56 percent strain. 


INTRODUCTION 

One important area of research m advancing the technology of aircraft gas 
turbine engines is to provide an accurate assessment of the life prediction procedures 
for blades and vanes of gas turbine engines . In order to further develop and eval- 
uate life prediction methods , simulated hardware components are tested in the labor- 
atory using carefully controlled conditions . Comparison of the experimentally meas- 
ured life to that which is analytically predicted is used as a means of evaluating life 
prediction theories . 

One experimental laboratory method for measuring thermal fatigue life is the 
cycling of wedge specimens in fluidized beds . Such tests can provide life data as 
well as transient temperature data that is obtained under carefully controlled condi- 
tions . Reference 1 contains a compilation of such data including a description of the 
facility and test procedure . References 2-9 contain incremental portions of such 
data relative to the evaluations described in this paper . 

Two fluidized beds were used for rapidly heating and cooling specimens in 
the form of prismatic bars with single-wedge constant cross-sections . These speci- 
mens failed by thermal fatigue cracking which is usually the predominant failure 
mode of aircraft engine first stage turbine blades and vanes . Thermal fatigue is de- 
fined as the cracking of a material caused by repeated temperature changes which in- 
duce cyclic stresses and strains . 

The turbine component life prediction methods currently being studied at Lewis 
are discussed in references 10-15. To predict specimen thermal fatigue life (defined 
as time-to-cracking) for the fluidized bed tests , it is necessary to know the stress/ 
strain-temperature-time history at the crack initiation site 
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The objective of this investigation was to analytically determine the elastic 
stress/strain-temperature-time history at the critical location for a single-edge 
wedge geometry specimen cycled in fluidized beds . This was performed as a 
joint program between NASA and the Air Force and utilized three-dimensional 
finite-element elastic analysis techniques . Lewis used the NASTRAN computer pro- 
gram while the Air Force Aero Propulsion Laboratory (AFAPL) used the IS03DQ 
program . 

Seven combinations consisting of three alloys and four different cycling con- 
ditions were analyzed. NASA TAZ-8A alloy was analyzed for the condition involving 
cycling in fluidized beds maintained at 316° and 1088° C (600° and 1990° F) , with 
the immersion time in each bed being 3 minutes. Alloys 316 stainless steel and 
A-286 were each analyzed for three different cyclic conditions. The heating bed for 
the three conditions was 740°, 850°, or 960° C (1364°, 1562°, or 1760° F) . The 
cooling bed was always 21° C (70° F) . Immersion time in both hot and cold beds was 
4 minutes . The cycling test conditions were those that resulted in thermal fatigue 
cracking in a reasonable number of cycles . The most severe testing condition was 
used for the high temperature alloy (NASA TAZ-8A) with less severe cycling condi- 
tions used for the other two intermediate temperature alloys . 

Due to symmetry , a discretized model of only a quarter of the wedge specimen 
was necessary for analysis . First a model with a fine mesh for a severe condition 
was analyzed using the NASTRAN program . Then a model with various coarse 
meshes for the same condition were analyzed using the IS03DQ computer program . A 
coarse mesh model for the IS03DQ analysis was selected which gave essentially the 
same results as using the fine mesh model with the NASTRAN analysis . The remain- 
ing combinations were then analyzed using the coarse mesh model and IS03DQ pro- 
gram . Such analyses provide the strainrange and stress/strain-temperature-time 
history so important for evaluation of life prediction theories such as strainrange 
partitioning . 

This work was conducted in the U.S. customary system of units. Conversion 
to International System of Units (SI) was made for reporting purposes only. 
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INPUT FOR COMPUTER PROGRAMS 

The alloys and test conditions for the seven combinations are given in table I . 
The necessary inputs to perform the analyses were: (1) the geometry of the single- 

edge wedge, (2) the elastic and physical material properties of the three alloys, and 
(3) a complete temperature distribution at various times throughout the cycle This 
section gives a detailed description of these inputs . 


Wedge Geometry 

The geometry for the single-edge wedge is shown m figure 1(a) The computer 
plots of the models used for analysis as well as a typical element are shown in fig- 
ure 1(b) for the NASTRAN program and m figure 1(c) for the IS03DQ program The 
model for both programs "squared-off" the leading-edge radius to a 1 02 mm (0 . 040 
in.) length. Otherwise the models duplicated the geometry of the wedge exactly. 
Detail discussion of the modeling is given m the section Description of Analyses . 


Alloy Properties 

The temperature independent and temperature dependent alloy properties used 
for the elastic analyses are given m tables II and III respectively . The properties re- 
quired for the analyses were Poisson's ratio, modulus of elasticity, and the mean 
coefficient of thermal expansion. The programs required a value for density to obtain 
results (zero mass elements were not permitted) although the results are independent 
of density . The properties for all alloys except the mean coefficient of thermal ex- 
pansion for NASA TAZ-8A alloy were obtained from references 16 and 17. The mean 
coefficient of thermal expansion for NASA TAZ-8A alloy was independently determined 
and is reported here for the first time. This, as well as all data in reference 16, was 
determined from the same heat used for fabricating the single-edge wedge test and 
calibration specimens . 
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Temperature Loading 

The transient temperature loading on the single-edge wedges was determined 
from thermocouple data . Calibration specimens of the three alloys were instru- 
mented chordwise at the midspan with five embedded thermocouples and cycled in 
the fluidized beds (schematically shown in figure 2) . The location of the thermo- 
couples at the wedge cross section is shown in figure 3 . The Inconel 600 sheathed 
thermocouples were mounted in grooves milled m the surface of the specimen and 
secured by a ceramic cement. The grooves were 0.56 mm (0.022 in.) wide and 
0.5 mm (0.02 in.) deep. Other details of the installation and procedure are given 
in reference 1 . The thermocouple outputs were cross-plotted to give temperatures 
of the midchord at the midspan at various time increments after immersion into the 
fluidized beds . These data are presented as figure 3 for the seven combinations 
analyzed . It was assumed that there was no temperature gradient through the thick- 
ness of the wedge . 

Another set of thermocouple data was taken with five thermocouples mounted 
along the leading edge over half the span . These data revealed a longitudinal 
(along the span of the blade) temperature gradient which varied with the different 
time increments . The maximum variation was about 16 percent greater at the ends 
of the wedge compared to the midspan and occurred after 30 seconds of heating . 
However , for any one time increment it was found that the ratio of the leading edge 
midspan temperature to that of any other span location was nominally the same for 
the seven investigated combinations . A least squares best fit parabola was deter- 
mined for each time increment and this is presented in table IV . This parabolic 
temperature variation along the span was assumed over the complete chord of the 
wedge . 

The temperatures at midspan were determined from the appropriate plot in 
figure 3. For locations other than midspan, the temperatures were determined by 
using the midspan temperature modified by the values given in table IV . There- 
fore, by use of figure 3 and table IV the temperature distribution at any point of the 
wedge was determined . 



DESCRIPTION OF ANALYSES 


Both computer programs used three-dimensional finite-element procedures 
to obtain an elastic analysis of the single-edge wedge geometry specimen . The 
NASTRAN program was used to obtain an analysis only for NASA TAZ-8A alloy 
for the time increment 15 seconds after immersion into the heating bed . The 
IS03DQ analysis , for the combination involving a 3 minute immersion time (fig- 
ure 3(a)) , was performed at all 17 heating and cooling time increments. These 
results showed that after 90 seconds immersion, the rate of change of stresses and 
strains was reduced to a point where coarser time increments might be used without 
loss of accuracy . Therefore , for the other six combinations which involved an addi- 
tional minute of immersion time, the analysis was performed at a total of only 16 
heating and cooling time increments - not using the 105, 135, 165, 195, and 225 
second time increments . 

The NASTRAN program was developed by NASA for elastic analysis of gen- 
eralized structures . The IS03DQ program was developed under contract by the 
Air Force for elastic analysis - specifically aircraft gas turbine blades , vanes , 
and disks . Documentation of the NASTRAN program includes a theoretical manual 
(ref. 18), a programmer's manual (ref. 19), a user's manual (ref. 20), and a demon- 
stration problem manual (ref. 21) . Documentation of the IS03DQ program includes 
a descriptive report (ref. 22) and a user's manual (ref. 23) . For general infor- 
mation on the programs the reader is referred to these manuals . Specific infor- 
mation on how the wedge was modeled and analyzed using these programs is pre- 
sented in the following sections 


NASTRAN Computer Program 

The model for the single-edge wedge was one-fourth of the structure as shown 
in figure 1(b) . There are reflective planes of symmetry at the midchord and mid- 
span for this structure . The nodal constraints on this model (using the axis nota- 
tion given in fig. 1(b)) are as follows 

1. No z-displacement for nodes on the midspan plane because of reflective 
symmetry . 

2. No -displacement for nodes on the midchord plane because of reflective 
symmetry . 
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3. Nox-djsplacemem for the two nodes at the trailing edge of the midspan 
plane to obtain a reference for displacements. 

A fine mesh was used in the NASTRAN analysis so that it might be used as the "base- 
line” for comparison. The model consisted of 464 nodes for the 192 CHEXA2 (hexa- 
hedral) elements. A typical element is shown in figure 1(b) . The discretization, 
including element and nodal identification, was done by hand - no mesh generator 
was used. The geometry was entered into the computer program by listing the co- 
ordinates of each node point from the origin shown in figure 1(b) Elements were 
selected so that the maximum aspect ratio for any element was always less than 2.5. 

The complete table of temperature dependent properties (modulus of elasticity 
and mean coefficient of thermal expansion) for NASA TAZ-8A alloy was entered with 
56° C (100° F) increments as given in table III. The program selected the value for 
these properties for each element by using the element temperature to linearly inter- 
polate within this table . Since temperatures were assigned to nodes rather than 
elements, the element temperature was determined by the program by averaging the 
eight nodal temperatures . Since NASTRAN does not have the capability to input 
temperatures by use of equations , all temperatures were first hand calculated (using 
fig. 3(a) and table IV) and then entered for each node point. 

The output selected from the NASTRAN program were the displacements , 
single point constraint forces , and stresses . The displacements and forces were 
given at the node points and the stresses were given at the element centroids. 
Stresses at the leading and trailing edges were obtained by extrapolation of plots 
through the centroids of the elements . 

This program was run using level 16.0 of NASTRAN on a Univac 1110 com- 
puter. The plot given in figure 1(b) was done on a Calcomp plotter using the 
NASTRAN plot subroutine . 


ISC jDQ Computer Program 

The model for the ISOSDQ analy~_ was similar to that used for the NASTRAN 
program. One-fourth of the single-edge wedge was used considering midchord and 
midspan planes of symmetry as shown in figure 1(c) . The nodal constraints were 
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identical to those used in the NASTRAN program so that a valid comparison could be 
made. The coarse mesh model selected consisted of 126 nodes for the 24 isopara- 
metric elements. A typical element is shown in figure 1(c) . The elements had mid- 
point nodes along the x-direction but not the y- and z-directions so that each ele- 
ment consisted of twelve nodes . The discretization , including element and nodal 
identification, was done using a mesh generator. This preprocessor (MESH 3) is 
part of the IS03DQ family of programs . This program required only the cross- 
section geometry of the wedge and some mesh parameters for the geometry input . The 
maximum aspect ratio for the elements was less than 17. 

Values for the two temperature dependent properties (modulus of elasticity 
and mean coefficient of thermal expansion) were entered into the program as segments 
of table III . This table gives the modulus of elasticity and mean coefficient of thermal 
expansion for each alloy at 56° C (100° F) temperature increments. Six of a program 
limit of eight values of modulus and thermal coefficient for six given temperature 
increments (table III) were put into the program . The program selected the value 
for the two temperature dependent properties for each node by using the nodal tem- 
perature to linearly interpolate within the table . 

The temperature loading was entered by means of a temperature table of eight 
chord temperatures at four different span locations . The program assigned a tem- 
perature to each node by weighted interpolation . Because temperatures were 
assigned to nodes rather than elements, a straight line gradient between adjacent 
nodes was assumed. 

The output selected from the IS03DQ program were the displacements , strains , 
and stresses . All of these values were determined at the node points . This set of 
data was put on tape for use by another program called PROUT3. The latter program, 
part of the IS03DQ family, allows the amount and format of the output to be varied 
without requiring the complete program to be rerun. Both the MESH3 (preprocessor) 
and PROUT3 programs have plot capability 

The IS03DQ family of programs were run on the Wnght-Patterson Air Force 
Base CDC 6600 computer. Plots (including figure 1(c)) were done using a Calcomp 
on-line plotter . 
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RESULTS AND DISCUSSION 

The results are presented and discussed in three parts. First, comparison 
of the NASTRAN and IS03DQ analyses for the check case are presented. Second, 
results for all seven combinations calculated by the IS03DQ program are presented 
at the critical location . The critical location was taken as that point on the blade 
which had the maximum longitudinal strainrange (algebraic difference between 
maximum and minimum longitudinal strain) throughout the complete heating and 
cooling cycle. This location was on the leading edge but not at midspan because 
of the longitudinal temperature gradient . Lastly , detailed computer plots for the 
seven combinations are presented at the times of both maximum and minimum 
longitudinal strain 


Comparison of NASTRAN and IS03DQ Analysis 

The comparison of the analyses of the single-edge wedge using the fine 
mesh model (fig. 1(b)) and NASTRAN with the coarse mesh model (fig. 1(c)) and 
IS03DQ is given in figure 4. The comparison shows very good agreement. Both 
analyses were independently performed for NASA TAZ-8A alloy after 15 seconds 
of fluidized bed heating . This alloy and time increment were selected as being 
approximately the most severe combination of all those studied to accentuate any 
differences between analyses. 

Figure 4(a) gives the normal x- , y- , and z-displacements along the leading 
and trailing edges . The greatest difference between the normal displacements as 
determined by the two methods is the x-displacement with the analyses for the other 
two displacements almost coinciding . However , comparing the net change between 
the leading and trailing edge x-displacements for the two analyses shows very 
similar results. 

Figure 4(b) gives the longitudinal stress along the midchord at about one- 
quarter span which was the critical location for this combination . These very good 
comparative results show that both the leading and trailing edges are in com- 
pression . This is due to the manner of testing in that the specimens were stacked 
so that they were heated and cooled from both the leading and trailing edges . A 
force balance of this cross-section showed that equilibrium requirements were 
satisfied . 
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This comparison confirmed that the IS03DQ program using a coarse mesh 
model was sufficiently accurate to obtain very good quantitative results . It also 
gave confidence in the use of this specialized blade and disk stress analysis 
program . 


Critical Locations 

Results for the seven analyzed combinations at the two critical locations 
(symmetrical about midspan) as a function of time after immersion into the 
fluidized beds are given m figure 5. This figure shows the temperature, and 
longitudinal strain and stress as a function of cycle time which occur at both criti- 
cal locations on the leading edge . 

In figure 5 , the temperature is the nodal temperature at the critical locations 
on the leading edge as determined from the temperature loading that was input to 
the IS03DQ program. The procedure used to determine this temperature is given 
in the section IS03DQ Computer Program . 

Both the longitudinal leading edge stress and strain show very steep gradients 
for about the first 10 seconds of immersion in both the heating or cooling beds . The 
results show that the leading edge goes into compression upon immersion into the 
heating bed. As the specimen reaches a steady-state condition, the stresses and 
strains approach zero. Upon immersion into the cooling bed, the leading edge 
goes into tension followed by a gradual drop-off to low stress and strain by the end 
of the cooling cycle . 

The maximum longitudinal strainrange for the combinations analyzed varied 
from 0.35 to 0.56 percent. NASA TAZ-8A alloy demonstrated the highest strain- 
range of 0.56 percent. For 316 stainless steel the strainrange varied from 0.35 to 
0.51 percent as the heating bed temperature increased. All three examined combin- 
ations of A-286 alloy exhibited strainranges of about 0.37 percent. 

Due to symmetry , the analysis showed two critical locations on the leading 
edge which were either 2.5 or 3.33 cm (1 .0 or 1 .3 in . ) away from midspan for 
each of the seven combinations evaluated . Preliminary experimental data in 
fluidized bed tests for some of the combinations indicate that cracks are initiated 
in this region . 
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The leading edge of the single-edge wedge is in a uniaxial state of stress . 
On the free surfaces, the normal x- and y-stresses (refer to fig. 1(c) for the 
axes convention) are zero. Therefore, the effective stress at the leading edge is 
equal in magnitude to the longitudinal z-stress . The x- and y-strains at the 
leading edge equal: 


where e = strain in x- , y- , or z-direction, and v = Poisson's ratio. 
By definition (ref. 24) effective strain is 




+ («> 


i ) 2 + (e c 


e l>' 


(2) 


where 1,2, and 3 refer to the principal directions . Since the shear strains are 
zero at the leading edge , the normal strains equal the principal strains . Substi- 
tuting eq. (1) in eq. (2) gives the effective strain at the leading edge as: 


„ - 2(1 + v) 

e eff 3 e z 


(3) 


Minimum and Maximum Longitudinal Strain 

Results for the seven analyzed combinations at the time increment of minimum 
and maximum longitudinal strain are shown in figures 6 and 7 respectively . The 
complete distribution of temperature and also normal , shear , and effective stresses 
and strains are shown over the complete midchord plane of the wedge . The nota- 
tion used is conventional elasticity notation with the axes convention as given in 
figure 1 . The assumption of constant temperature through the thickness of the 
wedge results in zero y-stress over the midchord . For this reason the y-stress 
plot is not presented . The minimum (largest compressive) longitudinal strain 
always occurred during heating and the maximum longitudinal strain always 
occurred during the cooling part of the cycle . These plots were made utilizing 
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the PR0UT3 program. These results show the reflective symmetry about midspan. 
These plots in addition to those in figure 5 will be used for further evaluation of 
various life prediction theories such as stramrange partitioning . 


SUMMARY OF RESULTS AND CONCLUSIONS 

The elastic stress analyses for a single-edge wedge geometry specimen cycled 
in fluidized beds were determined using three-dimensional finite-element techniques . 
The analyses were performed as a joint program of the NASA Lewis Research Center 
and the Air Force Aero Propulsion Laboratory (AFAPL) . The NASA TAZ-8A alloy 
was analyzed using the NASTRAN computer program with a fine mesh for a severe 
heating time increment. The AFAPL used the IS03DQ program with a coarse mesh 
model for this combination as well as all other combinations . A total of seven com- 
binations consisting of three alloys and four cycling conditions were considered . The 
alloys were NASA TAZ-8A, 316 stainless steel, and A-286. The cycling condition 
for NASA TAZ-8A alloy was alternate 3 minute immersions in cold and hot beds main- 
tained at 316° and 1088° C (600° and 1990° F) respectively. The three cycling con- 
ditions for both other alloys were alternate 4 minute immersions in a cooling bed 
always maintained at 21° C (70° F) and a heating bed maintained at one of three 
temperatures The temperatures were 740°, 850° , or 960° C (1364° , 1562° , or 
1760° F). Specific major results are as follows: 

1 . The analyses showed the leading edge of the single-edge wedge goes into 
compression when immersed into the heating bed followed by tension when immersed 
into the cooling bed . Steep stress and strain gradients occurred during the first 

10 seconds of immersion in either bed. For example, 0.35 percent strain was noted 
for NASA TAZ-8A alloy during the initial six seconds immersion in the heating bed. 

2 . The maximum longitudinal strainrange (algebraic difference between maxi- 
mum and minimum longitudinal strain) for the seven combinations considered varied 
from 0.35 to 0.56 percent. 

3 . The two locations of the maximum longitudinal strainrange at the leading 
edge of each wedge were either 2.5 or 3.3 cm (1 .0 or 1.3 in.) away from midspan 
for the seven combinations examined . Experimental test data for the combinations 
that have cracked indicate that the cracks initiated at these locations . 
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4. The comparison of the analyses using a fine mesh model (192 elements) 
and NASTRAN with a coarse mesh model (24 elements) and IS03DQ showed very 
good agreement for the single condition checked . 

5 . The results from this investigation can be used for further evaluation 
of various life prediction theories such as strainrange partitioning . 
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TABLE I. - ALLOYS AND CONDITIONS ANALYZED 


Alloy 

Fluidized bed temperatures 

Immersion time 
m each bed, 
sec 

Heating 

Cooling 

°C 

°F 

°C 

°F 

NASA TAZ-8A 

1088 

1930 

316 

600 

180 

316 stainless steel 

740 

1364 

21 

70 

240 


850 

1562 








960 

1760 







A- 286 

740 

1364 








850 

1562 








960 

1760 








TABLE II - TEMPERATURE INDEPENDENT 
ALLOY PROPERTIES 





















TABLE m - TEMPERATURE DEPENDENT ALLOY PROPERTIES 



NASA TAZ-8A 

316 stainless steel 


P 

l-286 


E a 

a 

b 

E a 

a 

b 

E a 

a b 

°c 

°F 

N/m 2 

psi 

m/m 

°C 

in. /in 
°F 

N/m 2 

psi 

m/m 

°C 

in /in 
°F 

N/m 2 

psi 

m/m 

°C 

in /m 
°F 

-17 8 

0 

208xl0 9 

30. 1x10 6 

10. 3xi0“ 6 

5 7x10" 6 

199xl0 9 

28. 8xl0 6 

15 3x10“ 6 

8 5x10“ 6 

201X10 9 

29 2x10® 

16 2xi0” 6 

9 0x10” 6 

37. 8 


206 

29.9 

10.6 

5.9 

197 

28 5 

15 5 

8 6 

199 

28 8 

16 4 

9 1 

93 0 


205 

29. 8 

11.0 

6 1 

194 

28 1 

15 7 

8 7 

195 

28 3 

16 6 

9 2 

149 


205 

29.7 

11.5 

6.4 

190 

27 6 

15 9 

8 85 

190 

27 5 

16 7 

9 3 

204 

400 

203 

29.5 

11.9 

6 6 

186 

27.0 

16 2 

9 0 

186 

27 0 

16 9 

9 4 

260 

500 

202 

29.3 

12. 1 

6.7 

181 

26.2 

16 6 

9 2 

183 

26 5 

17 0 

9 45 

316 

600 

201 

29. 1 

12. 1 

6 7 

177 

25.6 

16 9 

9.4 

178 

25 8 

17 2 

9 55 

371 

700 

199 

28.9 

12.2 

6.8 

172 

24 9 

17 3 

9 6 

174 

25 2 

17 4 

9 65 

427 

800 

198 

28.7 

12 4 

6 9 

166 

24 1 

17 6 

9. 8 

170 

24 6 

17 6 

9 75 

482 

900 

197 

28. 5 

12. 6 

7.0 

162 

23 5 

18 0 

10 0 

164 

23 8 

17 6 

9 8 

538 


194 

28.2 

12. 8 

7 1 

157 

22.7 

18 4 

10 2 

160 

23 2 

17 8 

9 9 

593 

1100 

192 

27.9 

12. 8 

7. 1 

152 

22 0 

18. 5 

10 25 

155 

22 5 

18 2 

10 1 

649 

1200 

190 

27. 5 

13 0 

7.2 

145 

21.0 

18. 5 

10.3 

152 

22 0 

18 4 

10 2 

704 

1300 

187 

27 1 

13. 1 

7 3 

139 

20. 2 

18.7 

10 4 

147 

21 3 

18.5 

10 25 

760 

1400 

183 

26. 5 

13.3 

7 4 

137 

19 8 

18 9 

10 5 

141 

20 5 

18 5 

10 3 

816 

1500 

178 

25. 8 

13.5 

7 5 

132 

19.2 

19 0 

10 55 

134 

19. 5 

18.7 

10 4 

871 

1600 

168 

24.3 

13.9 

7 7 

130 

18.9 

19. 1 

10. 6 

130 

18 8 

18 9 

10 5 

927 

1700 

146 

21 2 

14. 2 

7 9 

128 

18.6 

19 2 

10 65 

124 

18 0 

19 1 

10 6 

982 

1800 

139 

20.2 

14.6 

8. 1 

126 

18 3 

19 3 

10 7 

117 

17 0 

19 3 

10 7 

1038 

1900 

133 

19.3 

14.9 

8.3 

124 

18.0 

19 4 

10 75 

110 

15 9 

19 4 

10 8 

1093 


128 

18. 5 

15.3 

8.5 

123 

17 8 

19.4 

10. 8 

102 

14 8 

19. 6 

10 9 


a E modulus of elasticity. 


b a mean coefficient of thermal expansion from room temperature to indicated temperature. 
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TABLE IV. - TEMPERATURE VARIATION ALONG SPAN 


o 

[T = T (Az + Bz + C), where T _ is the temperature at any x, z 

L Ay 4 Ay UlO Xy Z 

coordinate (see fig. 1), T^ ms is the temperature at the x coordinate at 
midspah, and z is the span coordinate, all temperatures m °F (F = 

9/5 C + 32).] 


Time increment, 
sec 

Heating bed 

Cooling bed 

A 

B 

C 

A 

B 

C 

0 

-0. 00870 

0.0517 

0.9205 

-0. 00666 

0.03957 

0.9427 

3 

. 04401 

-.2614 

1.3891 

-.01775 

. 1055 

. 8447 

6 

. 03739 

-.2221 

1. 3290 

-.02384 

. 1416 

.7911 

9 

.03688 

-.2191 

1. 3372 

-.02548 

. 1514 

.7786 

12 

,03806 

-.2261 

1. 3344 

-.02731 

. 1622 

.7622 

15 


-.2195 

1. 3300 

-.02889 

. 1716 

.7480 

30 

sisa 

-. 1638 

1. 2504 

-. 03047 

. 1810 

.7338 

45 

.01769 

-.1051 

1. 1630 

-.03141 

. 1866 

.7224 

60 

. 01432 

-.08506 

1. 1324 

-.03442 

.2044 

.6905 

75 

. 01006 

-.05978 

1. 0934 

-.03265 

.1939 

.7093 

90 

. 00833 

-.04948 

1.0791 

-.02867 

. 1703 

.7440 

105 

. 00557 

-.03311 

1.0528 

-.02445 

. 1452 

.7843 

120 

. 00627 

-.03722 

1.0571 

-.02276 

. 1352 

.7981 

135 

. 00440 

-.02614 

1. 0415 

-.01876 

. 1142 

. 8323 

150 

.00371 

-.02205 

1. 0357 

-.01533 

. 09107 

. 8622 

165 

. 00297 

-.01762 

1.0285 

-.01278 

.07593 

. 8832 

180 

. 00262 

-.01553 

1.0243 

-.01212 

.07198 

.8876 

195 













210 













225 













240 
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Figure 2. - Schematic of fluidized bed test facility. 
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Figure 3 - Temperature of midchord at midspan at various times after immersion into the fluidized beds 
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Figure 5 -Continued 
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Figures -Continued 




Figures - Concluded 




r xy , 10 -4 m/m y xz , 10" 4 m/m y yz , 10' 4 m/m r max . 10 " 4 m/m 

(a) NASA TAZ-8A alloy with fluidized beds maintained at 1088° and 316° C 
(1990° and 600° F) after 9 seconds immersion in the heating bed 


Figure 6 - Temperature, stress, and strain distribution of midchord at time 
of minimum leading edge longitudinal strain (F - 9/5C + 32)11 ksi * 

6 895x10° N/m 2 ) 
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(c) 316 stainless steal alloy with the fluidized beds maintained at 850° and 21° C 
11562° and 70P F) after 9 seconds Immersion In the heating bed 


Figure 6, - Continued 
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id) 316 stainless steel alloy with the fluidized beds maintained at 960° and 21° C 
1 176(P and 7cP F) after 9 seconds immersion in the heating bed 

Figure 6 -Continued 
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(f) A-286 alloy with the fluidized beds maintained at 850^ and 21° C (1562° and 
7 (fi F) after 6 seconds immersion in the heating bed 

Figure 6 - Continued 
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(a) NASA TAZ-8A alloy with fluidized beds maintained at 1088° and 316° C 
(199CP and 600° F) after 9 seconds immersion in the cooling bed 

Figure 7 -Temperature, stress, and strain distribution of midchord at time of 
maximum leading edge longitudinal strain (F ■ 9/5C + 32)(1 ksi * 6 895X10 6 
N/m 2 ) 
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(b) 316 stainless steel alloy with the fluidized beds maintained at 740° and 21° C 
(1364° and 7(fi F) after 12 seconds immersion in the cooling bed 


Figure 7 -Continued 
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(c) 316 stainless steel alloy with the fluidized beds maintained at 850P and 21° C 
(1562° and 70° F) after 12 seconds immersion in the cooling bed. 



Figure 7 - Continued 










(d) 316 stainless steel alloy with the fluidized beds maintained at 96CP and 
21° C (1760° and 70° F) after 15 seconds immersion in the cooling bed 
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